Anyonic statistics with continuous variables 
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We describe a continuous-variable scheme for simulating the Kitaev lattice model and for de- 
tecting statistics of abelian anyons. The corresponding quantum optical implementation is solely 
based upon Gaussian resource states and Gaussian operations, hence allowing for a highly efficient 
creation, manipulation, and detection of anyons. This approach extends our understanding of the 
control and application of anyons and it leads to the possibility for experimental proof-of-principle 
demonstrations of anyonic statistics using continuous- variable systems. 



PACS numbers: 

I. INTRODUCTION 

It is well known from statistical physics that in three 
spatial dimensions, only two kinds of particles exist, 
bosons and fermions. A double exchange of any of these 
particles, topologically equivalent to braiding one around 
the other, will have no effect in three dimensions. In 
two spatial dimensions, however, braiding one particle 
around the other is no longer equivalent to doing noth- 
ing. Hence, the quantum mechanical wave function can 
acquire any phase through the braiding operation, lead- 
ing to the fractional statistics of exotic quasiparticles, so- 
called anyons Such anyons may be realized in model 
systems that have highly entangled ground states with 
topological order 0, y, 0| • The investigation of anyonic 
fractional statistics is not only of fundamental interest. It 
also led to the new field of topological quantum compu- 
tation, where the intrinsic robustness of the ground state 
Hamiltonian against local perturbations offers naturally 
given fault tolerance for quantum information processing 

EI 0,0. 

The concept of anyons first emerged in connection with 
the fractional quantum Hall effect 0], which occurs in a 
two dimensional electron gas at low temperatures being 
subject to a perpendicular magnetic field. The observa- 
tion of anyonic features in these systems requires cool- 
ing to very low temperatures where only the ground and 
a few excited states are populated. At the same time, 
anyons emerge as spatially localized many-body states 
that can be created and transported through the many- 
body system via local operations. Although some signa- 
tures have been observed in the fractional quantum Hall 
systems @ , a direct observation of fractional statistics as- 
sociated with anyon braiding has not been achieved there. 
Nonetheless, various theoretical proposals exist [l(|, in- 
cluding a spin-model scheme based upon laser manipu- 
lation of cold atoms in an optical lattice [ll|. However, 
currently, this scheme appears to be beyond current ex- 
perimental capabilities. 

A very simple model Hamiltonian from which one can 
obtain abelian anyons was proposed by Kitaev for spin- 



\ (qubit) systems This surface code model involves 
suitable combinations of four-body interactions. As such 
interactions are still hard to achieve on the level of the 
ground-state Hamiltonian, a more practical approach to 
simulating the surface code model would be based upon 
the creation and manipulation of multi-party entangled 
graph states [12]. Instead of cooling the system to its 
ground state, the ground and excited states are created 
dynamically from graph states. By focussing on the un- 
derlying multi-party entangled states, this approach pro- 
vides a testbed for demonstrating anyonic properties in- 
dependent of the presence of a Hamiltonian, in a simpler 
and more efficient way compared to the traditional ap- 
proaches based on complex solid-state systems (though 
in order to obtain the full robustness needed for fault- 
tolerant topological quantum computation, the back- 
ground Hamiltonian defining the protected code space 
appears to be crucial). 

The graph-state-based protocol was demonstrated ex- 
perimentally with six-photon (l3l | and four-photon graph 
states [l4j]. As reflected by these experiments, both the- 
oretical and experimental research of anyonic fractional 
statistics was mainly concentrated on two-level or spin- 
\ systems. Here, we shall consider anyonic statistics in 
continuous- variable (CV) systems [l5l |. i.e., for quantized 
harmonic oscillators representing, in particular, quan- 
tized optical modes. In order to apply the graph-state- 
based approach [l2j to CV systems, we shall consider 
CV cluster and graph states [l6j]. The corresponding 
optical Gaussian states can be efficiently and uncondi- 
tionally generated using off-line squeezed-state resources 
and linear optics [HI, [lTj . Various four- mode CV cluster 
states have been experimentally realized already [H, . 
By adding a non-Gaussian measurement to the toolbox 
of Gaussian operations, i.e., homodyne measurements 
within the model of measurement-based quantum com- 
putation (2(| , universal quantum computation is, in prin- 
ciple, possible using CV Gaussian cluster states [2lj . 

Here, within the CV setting, we will show that the 
anyonic ground state (or surface code state) can be cre- 
ated using Gaussian resource states (off-line squeezed 
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states) and linear optics. Generating quasiparticle ex- 
citations and implementing the braiding and fusion op- 
erations only requires simple single-mode Gaussian op- 
erations; also detection and verification of the quasipar- 
ticle statistics (observation of the fractional phase) can 
be achieved via Gaussian operations including homodyne 
detection. 

The plan of the paper is as follows. First, in SecHH we 
will briefly introduce the formalism and notations used in 
the paper for describing CV quantum logic and computa- 
tion. In the following, we discuss a generation scheme for 
the CV anyonic ground state based on CV graph states 
(Sec. IIIip and how to achieve the anyonic excitations in- 
cluding the braiding operations (Sec. IIV|) . In Sec. |Vj we 
will describe a detection scheme to verify the fractional 
phase of the anyons; in this measurement scheme, the 
global phase acquired by the excited state through the 
braiding operation is turned into a phase-space transla- 
tion resulting in different output stabilizer operators (cor- 
responding to different position- momentum linear combi- 
nations) depending on whether a braiding loop has been 
applied or not. Finally, before the conclusion in Sec. lVIIl 
we discuss in more detail potential optical implementa- 
tions of our CV-based anyon realization (Sec. IVI|) . 

II. CONTINUOUS- VARIABLE COMPUTATION 



Pauli group and the linear symplectic group of all one- 
mode and two- mode squeezing transformations [22l | . The 
CV gate for switching between the position and the mo- 
mentum basis is given by the Fourier transform operator 
F = exp[i(7r/4)(i 2 + p 2 )}, with F\x) x = \x) p . This gate 
corresponds to a generalization of the Hadamard gate 
for qubits. The phase gate P(rf) = exp[i(r]/2)x 2 ] is a CV 
squeezing operation, in analogy to a Clifford phase gate 
for qubits |23j| . The qubit controlled phase (sign) gate 
is generalized to a controlled-Z, Cz = exp(zii <g) X2). 
This gate provides the basic interaction for two modes 
1 and 2; at the same time, it describes a quadra- 
ture quantum nondemolition (QND) interaction. The 
set {X(s),F, P(i]),C z ;s,r) G M.} generates the Clifford 
group. Transformations within the Clifford group do not 
form a universal set of gates for CV quantum compu- 
tation. However, Clifford group transformations (corre- 
sponding to Gaussian transformations mapping Gaussian 
states onto Gaussian states) together with any higher- 
order nonlinear transformation acting on a single mode 
(corresponding to a non-Gaussian transformation with a 
Hamiltonian at least cubic in the mode operators) form 
a universal set of gates [22I l2~ij . As we will show in the 
following, all operations needed to demonstrate the any- 
onic fractional statistics in CV systems are entirely based 
upon Clifford group (Gaussian) transformations and ho- 
modyne measurements (Gaussian operations). 



In order to consider quantum logic over continuous 
variables, the qubit Pauli X and Z operators are gen- 
eralized to the Weyl-Heisenberg (WH) group [22j , which 
is the group of phase-space displacements. This is a Lie 
group with generators x = (a + at)/\/2 (quadrature- 
amplitude or position) and p — i(a— a))/ \/2 (quadrature- 
phase or momentum) representing, for instance, a single 
quantized mode (qumode) of the electromagnetic field. 
These operators satisfy the canonical commutation rela- 
tion [x,p] = i (with h = 1). In analogy to the qubit Pauli 
operators, the single- mode WH operators are defined as 
X(s) = cxp(— isp) and Z(t) — exp(iii;) with s,t € M. 
The WH operator X(s) is a position-translation opera- 
tor, which acts on the computational basis of position 
eigenstates {|x);x S R} as X(s)\x) = \x + s); Z(t) is a 
momentum-translation operator, which acts on the mo- 
mentum eigenstates as Z{t)\p) = \p + t). These operators 
are non-commutative and obey the identity 



X(s)Z(t) 



^(t)X( S ). 



(1) 



In the following, in analogy to the qubit case, we de- 
note the WH operators as Pauli operators. The Pauli 
operators for one mode can be used to construct a set of 
Pauli operators {Xi(si), Zi(ti); i — l,...,n} for n-mode 
systems. This set generates the Pauli group C\. The 
Clifford group Ci is the group of transformations that 
preserve the Pauli group C\ under conjugation; i.e., a 
gate U is an element of the Clifford group if URU f e Ci 
for every R G C\. The Clifford group Ci for continuous 
variables corresponds to the (semidirect) product of the 



III. CONTINUOUS- VARIABLE GRAPH AND 
SURFACE CODE STATES 

A graph quantum state is described by a mathemati- 
cal graph, i.e., a set of vertices connected by edges [25j . 
A vertex represents a physical system, e.g. a qubit (dis- 
crete two-level system in a 2-dimensional Hilbert space) 
or a CV qumode (system with a continuous spectrum 
in an infinite-dimensional Hilbert space). An edge be- 
tween two vertices represents the physical interaction be- 
tween the corresponding systems. The preparation pro- 
cedure of CV cluster states 16[ is analogous to that 
for qubit cluster states: first, prepare each mode (or 
graph vertex) in a phase-squeezed state, approximating a 
zero-momentum eigenstate (the analogue of the Pauli- X 
'+1' eigenstate); then apply a QND interaction [the gate 
Cz = exp(i£j (X) it)] to each pair of modes {j, k) linked 
by an edge in the graph. Note that all Cz gates com- 
mute. The resulting CV cluster /gra ph state (and, more 
generally, any CV cluster-type state [17|) satisfies, in the 
limit of infinite squeezing [16j . 



9a = (Pa 



E 

b£N a 



x b ) -> 0, Va S G, 



(2) 



where the modes a e G correspond to the vertices of 
the graph of n modes and the modes b £ N a are the 
nearest neighbors of mode a. More precisely, this re- 
lation describes a simultaneous zero-eigenstate of the n 
position- momentum linear combination operators g a . 
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FIG. 1: (Color online). Surface code model and anyonic 
braiding with continuous variables, a, The measure- 
ment pattern to prepare the anyonic ground state \tp) from 
a two-dimensional CV cluster state, where x and p denote 
a single-mode measurement in the position and momentum 
basis, respectively, F is the Fourier transform for all remain- 
ing modes. b,The CV surface code model consists of the star 
operators A s (£) and plaquette operators Bf(r)). c, Appli- 
cation of an X(t) position-translation on a single mode in 
the horizontal edge yields a pair of m(t)-type anyons placed 
at the neighboring plaquettes. d, Application of X(t) on a 
single mode in the vertical edge yields a pair of m(— t)-type 
anyons. e, Application of Z(s) on a single mode yields a pair 
of e(s)-type anyons at the neighboring stars. /, Application 
of X(t) on the modes along a string would create m-type 
anyons, where two m-type anyons are created on the same 
plaquettes. These two anyons annihilate each other accord- 
ing to the fusion rule m(t) x m(— t) = 1. Thus, the two single 
strings are glued together with two anyons at the endpoints. 
When a part of a string forms a loop, it cancels (dashed). 
g, An m-type anyon follows a closed loop around an e-type 
anyon. 



The stabilizer operators G a (£) = exp(— i£g a ) = 
X aiOUbEN a Z b(0 with £ e K for CV cluster states 
are analogous to the n independent stabilizers G a = 
X a Y\ b( - N Zb for qubit cluster states [25j. Starting from 
a 2D CV cluster state, single-mode measurements of ev- 
ery second mode in either the position or the momen- 
tum basis, with a subsequent Fourier transform F for 
all the remaining modes, yields a new graph state \ip). 
The measurement pattern is illustrated in Fig. la. An ad- 
measurement removes the measured mode from the clus- 
ter and breaks all connections between this mode and the 
rest of the cluster, including momentum-displacements 



of the nearest modes with the measured result. Con- 
versely, the effect of a p-measurement is to remove the 
measured mode, but to connect the neighboring modes 
at the same time. These measurements transform the 
set of correlations in terms of position-momentum lin- 
ear combination operators from {{p a — YlbeN ®b) ~ * 0} 
to {a s = (p s> i + p S:2 + p s ,3 + Psa) -> 0, bf = - 
%f,2 + £/,3 — £f,i) ~ * ^} f° r the remaining modes, where 
s and / label stars and plaquettes (faces), respectively, 
and the indices 1,...,4 of the position and momentum op- 
erators denote those modes located at a common star or 
at the boundary of a common plaquette. The stabiliz- 
ers {A s (£) = exp(-i£a s ) = U jestar ^X s> j(^),Bf(ri) = 

exp(-W?6/) = Ujeboundary(f)Zf,j((- i y , n)}, witn fj V e 

M., for the new CV state are exactly analogous to 
those of the first Kitaev model Q for a two-dimensional 
spin lattice. These stabilizers, the star operators A s (£) 
and plaquette operators Bf{q), all commute. So the new 
CV state \tf)) corresponds to the anyonic ground state 
with MZM = \4>) and Bf(rj)\f) = |ifr>, VC,77 € E, 
for all stars s and plaquettes /, in the limit of infinite 
squeezing. After defining the CV anyonic ground state 
and describing how to create an approximate version of 
this state (asymptotically perfect in the limit of infinite 
squeezing) from CV cluster states, let us now consider 
anyonic excitations from the ground state and the basic 
braiding operations. 



IV. ANYON CREATION AND BRAIDING 

Starting from the ground state one can excite pairs of 
anyons connected by a string using single-mode opera- 
tions. Unlike the Pauli group for qubits, the WH group 
is a continuous group, and therefore quasiparticle exci- 
tation is continuous. More specifically, by applying the 
position-translation operator X(t) to some mode of the 
lattice, a pair of so-called m-type anyons in the state 
\m((-l) d t)) = X(t)\ip) (d G {1, 2}, d = 1 means that the 
relevant mode lies on the vertical edges, whereas d = 2 
refers to the horizontal edges) is created on the two neigh- 
boring plaquettes (Fig.lc,d). An e-pair of anyons, given 
by \e(s)) = Z(s)\i/j), is obtained on two neighboring stars 
by a momentum-translation operation Z(s) (Fig.le). 

Similar to the qubit case, the excited states differ 
from the ground states in their stabilizer operators. In 
the CV case, this means that some of the position- 
momentum linear combination operators a s and b f [from 
the ground-state stabilizers A s (£) = exp(— i^a s ) and 
Bf(rj) = exp(— irjbf)] no longer satisfy the ground-state 
conditions d s = and bf — (in the limit of infinite 
squeezing), but rather become a s = t for Zj 5star ^(t)\ip) 

and b f = ±s for X jeboundary ^(s)\ip). In other words, 
the anyonic excitations become manifest in "violations" 
of the stabilizer ( "nullifier" ) conditions for a s in the 
presence of e-particles and for bf in the presence of m- 
particles. 
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The fusion rules [e(s) x e(t) — e(s + t), m(s) x m(t) = 
m(s + t), e(0) = m(0) = 1, 1 x e(s) = e(s), 1 x m(s) = 
m(s)] describe the outcome from combining two anyons. 
Thereby, excited states can be obtained from the CV 
ground state \ip) by applying open string operators, 
which create quasiparticles at their endpoints (Fig. If). 
The closed X (Z) strings surround the corresponding 
stars (plaquettes) , which are a product of the surrounded 
star (plaquette) operators (Fig. If). Thus, for any closed 
X (Z) string, the system will remain in the ground state. 

The anyonic character of the excited states is now re- 
vealed through a non-trivial phase factor acquired by the 
wavefunction of the lattice system after braiding anyons, 
i.e., after moving m around e (Fig.lg) or vice versa. Con- 
sider the initial state l^mi) = Z\{s)\ijj) = |e(s)). If an 
anyon of type m is assumed to be at a neighboring plaque- 
tte, it can be moved around e along the path generated 
by successive applications of X(t) on the four modes of 
the star. The final state is 

\V fin ) = Xi(t)X 2 (t)X a (t)* 4 (t)|tf <n< ) 

= e- lst Z x (a) [X ± (t)X 2 (t)X 3 {t)X 4 (t)\if>)] 

= e- is *|^ m ), (3) 

using Eq.JI]). The extra factor exp(— ist) is the topo- 
logical phase factor, which reveals the presence of the 
enclosed anyons (it would not be obtained, if the ini- 
tial state had been the unexcited state, l^mi) — \tl>))- 
Note that here in the CV setting, the quantum state can 
acquire any phase, in contrast to the discrete value for 
anyons in qubit systems with the analogous braiding op- 
eration (where XZ = —ZX for qubits). 

The corresponding phase factors, changing the 
phase of the wavefunction, constitute an abelian (com- 
mutative) representation of the braid group. It is 
important to understand that, in general, a relation 
as described by Eq.([3]) can be obtained for different 
paths of the loop. For instance, the loop 1-2-5-6-7-4 
of m-type anyons (Fig.lg), which corresponds to the 
application of X 1 (t)X2(t)X 5 (-t)X 6 (-t)X 7 (-t)X 4 (t) = 
[Xti^^Xs^X^Xsi-t^i-^Xei-^Xri-t)}, 
involving two star operators, gives the same result for the 
topological phase as in Eq.([3]). In fact, this equivalence 
reflects the topological character of the extra phase 
factor and the potential robustness of the corresponding 
state transformation, as it is exploited in topological 
approaches to fault-tolerant quantum computation. 

V. ANYON DETECTION 

In order to detect the most significant feature of 
anyons, their non-trivial statistical phase obtainable 
through braiding, we may employ an "interference mea- 
surement" in phase space, which enables us to detect the 
otherwise invisible global phase factor of Eq.([3]). The 
idea is to turn the global phase acquired by the excited 
state into a displacement in phase space that would not 
occur without braiding operation. 



First, we generate a superposition of the ground state 
\tp) and the continuously-excited state of the anyon e 
via the squeezing operation P(s) = exp[i(s/2)x 2 ] (the 
CV quadratic phase gate) instead of the momentum- 
translation operation Z(s) — exp(isi). Then we perform 
a closed loop of anyon m around e. Finally, application of 
the inverse squeezing operation P(— s) = exp[— i(s/2)x 2 ] 
makes the phase difference visible, leading to the state 
e -ii -s/ 2 1 e (si) } according to 

e -i{s/2)x 2 e ~itp e i{s/2)x 2 _ ^-it 2 s /2 e istx & ~itp U\ 

This compares to the state e -*(»/ 2 ) a:2 e *( s / 2 ) x2 — |^,) 
without a closed loop of anyon m around e. In other 
words, the effect of the braiding operation can be mea- 
sured via the different stabilizer operators or, more pre- 
cisely, the different position-momentum linear combina- 
tion operators of the output states: for the relevant star 
operator so, we have a So = st with braiding and a So = 
without braiding. 

Note that this verification scheme is only complete, 
provided that, in addition to detecting the change of the 
stabilizer (nullifier) operator a So , a sufficient set of extra 
milliners a s and bf is measured in order to confirm the full 
inseparability of the output states (the ground state and 
the momentum-displaced ground state) [26j]. Otherwise, 
even a fully separable product state of zero-momentum 
eigenstates used as "ground state" could produce the 
same braiding statistics, as it would satisfy a So — with- 
out braiding and a So = st with a "closed braiding loop" . 
Such a trivial ground state, however, would obviously 
not simulate Kitaev's surface code model; in particular, 
apart from not satisfying the complete set of ground-state 
correlations, it would not allow for a verification of ex- 
cited TO-particles or the "conjugate" braiding operation 
of moving e around to. 

VI. OPTICAL IMPLEMENTATION 

In a real experiment, because of finite squeezing, it 
is impossible to create perfect correlations in terms of 
position-momentum linear combination operators (the 
stabilizers or, better, milliners) for the anyonic ground 
state l^)- Nonetheless, we may still implement and de- 
tect anyonic fractional statistics by interpreting the pres- 
ence of small fluctuations around the ground-state stabi- 
lizer conditions due to finite squeezing, without the ex- 
istence of any first-moment phase-space translations, as 
the ground state (Fig. 2d); sufficiently large first-moment 
position (or momentum) translations, still maintaining 
the second-moment multi-mode position-momentum cor- 
relations, would then correspond to the excited states. 

For verification, the amplitude (phase) quadrature, 
i.e., the position (momentum) of each mode is measured 
via homodyne detection and every detector output is 
sampled within a time interval to yield a number of mea- 
sured quadrature values. Using the correlation diagram 
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FIG. 2: (Color online). Surface code state with four 
modes used for demonstration of anyonic statistics. 

a, The anyonic ground state \iP)a, forming one star and four 
plaquettes, is given by a four-mode GHZ state with the modes 
labelled 1 to 4. b, The four-mode cluster state which after 
Fourier transforms on modes 1,3, and 4, is equivalent to the 
anyonic ground state \ip)i- c, The quantum circuit for cre- 
ation, braiding, and annihilation of the anyons, and for the 
final state detection, d, The correlation diagram of the mo- 
mentum linear combination operator pi + f>2 + P3 + Pa with 
(red) and without a closed loop of anyon m around e (black) 
in the case of finite squeezing. 



[27| containing a sufficient number of measured quadra- 
ture values, one can determine first and second moments 
of the position-momentum linear combination operators 
for the anyonic ground and excited states. Thus, the fully 
inseparable ground states can be verified and discrimi- 
nated from the (for large squee zing near-orthogonal) ex- 
cited (displaced ground) states [26[ ; at the same time en- 
abling one to confirm the braiding-dependent state tran- 
sition according to Eq.((4])(see Fig. 2d). For finite, even 
small squeezing values, it is still possible to discrimi- 
nate the ground and excited states, provided the ampli- 
tude (or phase) translations are sufficiently large. Higher 



squeezing values improve the sensitivity for detecting the 
amplitude (or phase) translations of the excited states. 
Note that even for zero-squeezing resources, such a state 
discrimination can be almost perfect (for large ampli- 
tudes/phases); in this case, more subtle is the verifica- 
tion of the full inseparability [26] of both the ground and 
excited states, as the multi-party entanglement in the ap- 
proximate CV surface-code model would decrease in this 
case (but it would not entirely vanish, because of the 
intrinsic squeezing of the Cz gates for CV graph-state 
generation) . 

In the Kitaev lattice model, the quasiparticles are par- 
ticularly well localized: an e-particle is on a single star 
and an m-particle is on a single plaquette. Based on 
this perfect localization, the proposal for implementing 
the anyonic model via the simple surface code has been 
demonstrated experimentally with four-photon [l4j and 
six-photon fl3j graph states. The four-photon experi- 
ment relies upon only four photonic qubits entangled in 
a four-party GHZ state. This simple system is sufficient 
for demonstrating a small Kitaev lattice (with one star 
and four plaquettes) and a non-trivial braiding operation, 
because the remaining eight qubits of the total 12-qubit 
system factor out from the four-qubit GHZ state in a fully 
separable product state [14j. In the CV regime, we pro- 
pose to pursue an analogous strategy. In fact, four- mode 
entangled CV GHZ and cluster states have been demon- 
strated already experimentally [HI, . In a small-scale 
surface-code experiment, it is crucial to verify the full in- 
separability of the GHZ-type ground state as well as that 
of the excited state, i.e., a displaced GHZ state. These 
simple systems may then also be used to implement and 
detect minimal braiding loops (Fig. 2a,b,c) with four- 
mode entangled GHZ states. The complete verification 
requires detecting the position-momentum linear combi- 
nations corresponding to the stabilizer set 

X 1 (t)X 2 (t)X 3 (t)X 4 {t), 

Z x {t)Z^{t), Z 2 (t)Z^(t), Z 3 (t)Z^(t), (5) 

including first and second moments [26| . 

However, for the simple four-mode GHZ states, there is 
no way to show the robustness of the topological braid- 
ing operation when the anyons are moved along differ- 
ent paths. In order to demonstrate such robustness of 
the anyon operations against different braiding paths, 
more modes would be needed, for instance, nine-mode 
states in analogy to the nine-qubit states of Ref.[l2j. 
The corresponding CV nine-mode ground state \ip)g is 
locally equivalent to a CV nine-mode graph state, which 
can be efficiently generated via nine off-line sing le-mode 
squeezed states and arrays of beam splitters [T7| instead 
of the projective measurements on a prepared 2D CV 
cluster state as described previously [28|]. In these larger 
lattices, there are several different loops for the quasi- 
particles. The anyonic statistics, however, would only 
depend on the topological character of the loop. 
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VII. CONCLUSION 

In summary, we have proposed a protocol to create, 
manipulate, and detect anyonic states, including demon- 
stration of their fusion rules and their fractional statis- 
tics, in the regime of continuous quantum variables. 
In our continuous-variable approach to Kitaev's surface 
code model, the main features of anyons can be obtained 
via simple Gaussian operations starting from the creation 
of suitable continuous-variable graph states. With re- 
gard to a potential optical realization, this means that 
all ingredients for a full verification can be implemented 
in a very efficient way: graph-state generation only re- 
quires off-line squeezing and linear optics; the surface- 
code ground states are then obtained either through 
homodyne measurements or via local phase rotations; 
anyon creation and braiding-dependent state transitions 
can be also detected through Gaussian operations, i.e., 
online s que ezing transformations (which could be shifted 
off-line [30] ) and homodyne detections (including first 
and second moments of quadrature combinations in order 
to verify the fully inseparable ground and excited states 
and to discriminate them against each other). 



This work is a first step towards further theoretical 
investigations of topological quantum computation over 
continuous variables and experimental proof-of-principle 
demonstrations of anyonic, topological features via rela- 
tively simple, stable optical systems, as opposed to the 
more traditional approaches within condensed matter 
physics. 
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